Electron spin relaxation in graphene: the role of the substrate 
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Theory of the electron spin relaxation in graphene on the SiC>2 substrate is developed. Charged 
impurities and polar optical surface phonons in the substrate induce an effective random Bychkov- 
Rashba-like spin-orbit coupling field which leads to spin relaxation by the D'yakonov-Perel' mecha- 
nism. Analytical estimates and Monte Carlo simulations show that the corresponding spin relaxation 
times are between micro- to milliseconds, being only weakly temperature dependent, ft is also argued 
that the presence of adatoms on graphene can lead to spin lifetimes shorter than nanoseconds. 

PACS numbers: 72.25.Rb, 73.61.Wp, 73.50.Bk 



Since the experimental realization of graphene, a sin- 
gle stable 2D-monolayer of carbon atoms arranged in a 
honeycomb lattice, considerable research has been done 
to enlighten its peculiar electronic transport properties 
originating from the Dirac-like band structure at the K 
and K' points in the momentum space [l| . Long spin re- 
laxation times and phase coherence lengths in graphene 
are expected based on the weak atomic spin-orbit cou- 
pling in carbon [Z = 6). However, recent spin injection 
measurements based on a non-local spin valve geometry 
0, 0, 0] revealed surprisingly short spin relaxation times 
of only about 100-200 ps, being only weakly dependent 
on the charge density and temperature. These results 
appear puzzling, although the low mobilities of the sam- 
ples (about 2000 cm 2 /Vs) suggest that the measured spin 
relaxation times are likely due to extrinsic effects . 

Very recent experiments on the charge transport in 
graphene affirmed the importance of the underlying sub- 
strate [1, @, 0| ■ At low temperatures the transport prop- 
erties have been shown to be dominated by scattering 
from the charged impurities residing in the substrate 
@, Q . The conductivity of graphene placed on a SiC"2 
substrate starts to decrease above 200 K. The observed 
temperature and density dependence of the resistivity are 
most likely explained by remote phonon scattering due to 
occurrence of polar optical surface modes in the substrate 

BEE!. 

These findings naturally raise the question if (i) 
charged impurities and (ii) remote surface phonons are 
also relevant for the spin relaxation in graphene. As 
argued here both mechanisms provide a temperature- 
dependent, random spin-orbit coupling field, which limits 
the spin relaxation via the D'yakonov-PereP (DP) mech- 



anism [II El ELI- The calculated spin relaxation times 
are micro to milliseconds. In addition, we give estimates 
for the spin relaxation times due to the possible presence 
of adatoms on graphene. For reasonable adatom densities 
the spin lifetimes can be lower than nanoseconds. 

Several other mechanisms have already been investi- 
gated theoretically, such as the spin relaxation due to the 
corrugations (ripples) of graphene and due to exchange 
interaction with local magnetic moments [HI], or spin- 
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FIG. 1: (Color online) (a) The Dirac cones when spin-orbit 
coupling is included. The arrows indicate the spin vectors 
rimv as defined in the text, (b) Effective magnetic field direc- 
tions (Bychkov-Rashba field) along the Fermi circle of elec- 
trons at the Jf-point (the field is the same at the 7f '-point), 
(c) Graphene layer on the top of a Si02 substrate with charged 
impurities which induce an electric field component E z per- 
pendicular to the plane breaking the inversion symmetry of 
graphene. (d) Illustration of the spin relaxation in a spatially 
random potential due to the charged carriers. In the Monte- 
Carlo simulations the spin dynamics is sampled over random 
trajectories with different initial momenta. 



orbit coupling mediated relaxation based on boundary 
scattering, heavy impurities, and effective gauge fields 
due to topological disorder [16| . 

Near the K and K' points the carrier dynamics can be 
described by an effective low energy Hamiltonian 17| of 
the form Hq = hvf(T z a x k x + <J y k y ). Here, Vf — 10 6 m/s 
denotes the Fermi velocity, k is the wave vector with re- 
spect to K(K'), and r and a are the Pauli matrices with 
t z = ±1 describing the states at the K and K '-points, 
respectively, and a z = ±1 describing the states on the A 
and B sublattice of the honeycomb lattice. The inclusion 
of the microscopic spin-orbit (SO) interaction results in 
an additional term in the effective low-energy Hamilto- 
nian, Hi = —\i + \it z o~ z s z as shown either by group the- 
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oretical arguments [l8[ or by second order perturbation 
theory of a microscopic tight-binding model ljj 2(3, 21 1 . 
Here, the real spin is represented by the s z Pauli ma- 
trix and Xj denotes the intrinsic SO-constant. The in- 
trinsic SOC opens a gap Aj = 2Aj at the Dirac point, 
making graphene theoretically a spin Hall insulator [lj| . 
Recent first-principles calculations give Aj = 0.024 meV 
[22l |. large enough to influence electronics of intrinsic or 
weakly charged graphene only somewhat below 1 K. 

If an electric field is applied perpendicular to the 
graphene plane, the inversion symmetry is lifted and 
group theory allows for an additional Bychkov-Rashba 
(BR) term of the form Hrr = Xbr(t z o~ x s % 



[la , Il9l . fell |23|. From first-principles calculations [25 
a linear relationship between the BR-constant and the 
electric field, Abr(i") = CbrE z (t), is found, with (br = 
0.005 mcV/(V/nm). The proper knowledge of Abr is of 
great importance for our quantitative results below, since 
in the DP-mechanism the spin relaxation rate depends 
quadratically on Abr. 

The resulting effective 8x8 Hamiltonian H c r = Ho + 
Hi + Hbr is easily diagonalized yielding the same eigen- 
values at the K and K' points, 
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v\i) + myje 2 + (A B r - vXi) 2 (1) 



with e = hvf\k\. v = ±1, and the band index m = 1 for 
electrons (e) and m = — 1 for holes (h), respectively. We 
define spin vectors vH mv = s^/lsJ^J as normalized ex- 
pectation values of the spin operator s^ nv = (ip\S\tp) with 
respect to the eigenstates \ip) = |r, k, m, v) (t = K,K') 
of the total Hamiltonian iJ c ff- The vectors are in-plane 
and result in rig + = nj l _ = (— sirup, cos tp, 0) and nl_ = 
nj l+ = (sin tp, — cos tp, 0) with tp denoting the polar angle 
of the wave vector k. In the case of e S> Ah + Aj, i.e., 
if the Fermi energy is much greater than w 0.03 meV (a 
condition usually fulfilled in gated or doped graphene) 
the electron and hole motion can be decoupled. By suc- 
cessive unitary rotation of H e g first into the eigenbasis 
of Hq and then into the spin basis with respect to the di- 
rection n = (— smtp, cos tp, 0) an effective BR-typc 2x2 
Hamiltonian can be obtained for both holes and elec- 
trons, 



H e B = m(s — Xi) + mABRn(k) ■ s 



(2) 



with s denoting the Pauli spin matrices. H e g is the same 
for K and K' , as guaranteed by time reversal symmetry. 

Comparison with the original BR-Hamiltonian in semi- 
conductor heterostructures of the form iJk = M~2(k) • s/2 
shows that SOC coupling in graphene effectively acts on 
the electrons spin as a in-plane magnetic field of con- 
stant amplitude but k-dependent direction, as illustrated 
in Fig. [Ha) an d (b). In this effective field the spin 
precesses with a frequenc y o f Q = 2Abr/S. As shown 
by D'yakonov and Perel' [l2l . [l3| random scattering in- 
duces motional narrowing of this spin precession causing 
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FIG. 2: (Color online) Calculated spin relaxation time r s as 
a function of the Fermi energy Ef, taking into account only 
impurity scattering, for two different impurity densities in the 
substrate at T = K. The symbols indicate MC-simulation 
results and the solid lines are analytic fits of the form 1 /t s = 
T imp (Ef)Q,gg with O e g = 3.3x 10 9 s _1 (for squares) and fl e g = 



1.1 x 10 8 s" 1 (for triangles) 



spin relaxation. The spin relaxation rates for the DP- 
mechanism for the a-th spin component generally result 
in l/r s , a = T *((^k) ~ (^a)) w i tn T * denoting the corre- 
lation time of the random spin-orbit field and (...) indi- 
cates averaging over the Fermi surface. Due to the po- 
lar angle dependence, in graphene the correlation time 
exactly coincides with the momentum relaxation time 
t* = T p [l3, 14 1 . Hence, for graphene the spin relaxation 



2r s 



time results in 1/t s _ z = t p (2Abr/S-) 2 and t. 

First, we investigate spin relaxation due to charged 
impurities residing in the substrate as schematically il- 
lustrated in Fig. []Jc). Impurity scattering is a dominant 
scattering mechanism governing the transport properties 
of graphene Due to the fluctuations of the impu- 

rity concentration a random unscreened electric field per- 
pendicular to the graphene plane, and hence a spatially 
random BR-field Abr(i") = (brE z (t) arises. As shown by 
Sherman 24j in the case of semiconductor quantum wells, 
the randomness of the BR-field in the real space already 
causes spin relaxation even without any scattering in the 
k-space. The correlation length of the random BR-field 
is on the scale of the distance d of the impurity layer from 
the graphene sheet j24[. Therefore, the spin relaxation 
time for a ballistically moving electron can be estimated 
as follows: if the electron passes through a domain of 
the lateral size of the correlation length of the BR-field, 
the spin precesses by 5tp = flQRd/vf. At some time t 
the electron has passed through t/(d/vf) different do- 
mains and in the picture of a random walk it follows that 
Stpy/t/ (d/vj). The spin is relaxed if ((A^) 2 ) w 1 



yielding the condition l/r s 
a semiclassical picture for the orbital motion r(t) of the 
electron, the spin experiences a random BR-field both in 



4/h 2 (X BR )d/vf. Hence, in 
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the real space (Sherman mechanism) and in the recipro- 
cal space due to momentum scattering (DP-mcchanism). 

We numerically calculate the spin relaxation time by 
performing Monte Carlo (MC)-simulations for the spin 
dynamics. For this purpose we use a random but 
quenched impurity distribution of a given density and 
sample over the random particles trajectories starting 
with different initial momenta, as illustrated in Fig.[TJd). 
The trajectories are generated according to the scatter- 
ing probability of the screened impurity potentials in the 
graphene sheet calculated in the random phase approxi- 
mation following Ref. @. 

Along any given semiclassical trajectory [r(t), k(f)] the 
spin dynamics can be described by the Bloch equation 



ds 

dt 



n BR [r(t)](n[k(t)] xs). 



(3) 



The spin relaxation time is then calculated by averaging 
over the asymptotics of all trajectories, since for times 
much greater than the mean free time t r m f p the spin 
components relax as s a (t) = s Q (0) cxp(— t/r St a) [Ill- 
Figure [2] shows the calculated spin relaxation time as 
a function of the Fermi energy Ef for a dirty SiG"2 sub- 
strate, n imp = 4 x 10 12 cm~ 2 , and for a cleaner sample, 
"imp = 4 x 10 11 cm , taking into account only impu- 
rity scattering. For all simulations we use the ab-initio 
BR-parameter (br = 0.005 meV/(V/nm) and an effec- 
tive impurity distance of d = 0.4 nm from the graphene 
layer p, |9j . The symbols refer to the MC-simulation re- 
sults and the solid lines indicate analytic fits of the form 
1/t s = T imp (E /)^cff with Ti m p denoting the momentum 
relaxation time due to impurities. Since the cross section 
of the screened long-ranged Coulomb potential is propor- 
tional to the Fermi wavelength A/ ~ kj 1 the mo- 
mentum relaxation time increases with increasing Fermi 
energy yielding a decreasing spin relaxation time, as il- 
lustrated in Fig. [2] 

The second important spin relaxation mechanism in- 
duced by the SiC>2 substrate is due to polar optical sur- 
face phonons. In the case of SiC"2 there are two dominant 
surface phonons with energies of tkoi^ = 59 meV and 



(2) 



155 meV, respectively, which provide a temper- 



ature dependent electric field variance given by [2E 



{El i )(T)=p, 



*JP (l + 2n«) 



47T£n 



4d 3 



(4) 



with £o denoting the dielectric constant and standing 
for the Bose-Einstein occupation factors of the phonon 
mode i. The individual strengths of these remote phonon 
scattering modes are given by the parameters fli = 0.025 
and (i 2 = 0.062, which fulfill the relation (3 = ^< A = 
(e s - Eoo)/{e s + 1)(£«> + 1) with (3 giving a measure 
of the total polarizibility of the dielectric interface Q 
and e s and Eoo denoting the static and high-frequency 



5 4 
2 



■1/t, 

■1/t. 

i" 

1/t 
1/t 



tot 



imp 
sph 



°0 



100 



200 



300 



400 



FIG. 3: (Color online) Calculated inverse momentum relax- 
ation times 1/t as function of temperature T for impurity 
(imp) scattering (with n lrap — 4x 10 11 cm -2 ), surface phonon 
(sph), and acoustic phonon (ac) scattering at Ef = 100 meV. 



dielectric constant, respectively Due to the random- 
ness of the electrons' motion the spin experiences an 
effective electric field and, hence, a random BR-field. 
The effective spectral correlation function of the phonon 
field (E z {t)E z (t')) will include an exponential decay with 
the momentum relaxation time r m yielding a Lorentzian 
rcnormalization factor E 2 ai = (£ 2 ^/[l + (a>^V m ) 2 ] [l4|. 
If uj s T m 3> 1 (as for graphene on Si0 2 ), the effective 



electric field can be found by qualitative arguments: for 
long-wave phonons the spin precesses by Sip = fiBR T ph 
in the characteristic time r p h = l/u) s . The momen- 
tum scattering leads to a random walk with typical 
step times of r m . The spin is relaxed if the variance 
(Sip) 2 = (t/T m )Sip reaches one, yielding for the spin re- 
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giving 



laxation time 1/t s = r2g R r m /(a;phT m ^ 
an effective field of E 2 a = E 2 / (uj p hT m j 2 ■ 

For the MC-simulations we took into account momen- 
tum scattering due to charged impurities in SiC>2, op- 
tical surface phonons and acoustic phonons of the 
graphene sheet [ll|. The resulting total momentum re- 
laxation rate l/r to t = 1/T imp + l/T sph (T) + l/r ac (T) is 
illustrated in Fig. [3l showing that the impurity scatter- 
ing remains dominant up to room temperature but with 
an exponentially increasing contribution coining from the 
surface phonons and a linearly growing contribution due 
to acoustic phonon scattering. The random BR-field is 
calculated from the electric field originating from the im- 
purities and the polar surface phonons. The temperature 
dependence of the spin relaxation time for a fixed Fermi 
energy of Ef = 100 meV for different impurity densi- 
ties is shown in Figure [H where the solid lines indicate 
again fits of the form 1/t s = T tot (T)[2(BR,E e ff / H] 2 . The 
MC-simulations reveal that the spin relaxation time is 
almost temperature- independent. This is caused by the 
nearly perfect counterbalancing of the increasing electric 
field and the decreasing momentum relaxation time with 
temperature. As for the relaxation of the momentum, 
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FIG. 4: (Color online) Calculated spin relaxation time t s 
versus temperature T taking into account impurity, surface 
phonon and acoustic phonon scattering at Ej — 100 meV. 
The symbols refer to MC-data and the solid lines are fits of 
the form l/r s = T tot (T)[2C, B nE ea / h] 2 with E cS = 0.21 V/nm 
(for squares) and E e g = 0.007 V/nm (for triangles). 



impurities dominate the spin relaxation compared to the 
mechanism of optical surface phonons, which causes a 
decrease of r s by about 10 — 20%. 

Can we relate our results to the experimental findings 
of t s of 100-200 ps 0? Even considering the uncertain- 
ties in d or in the charge density in the substrate, such 
small values for r s can be hardly explained by the sub- 
strate effects. Indeed, the measured samples have short 
mean free times of about T m f p w 50 fs 0] , which suggest 
a high impurity density of about n\ mp = 2 — 4 x 10 12 
cm -2 0. However, the times r s « 100 ps would require 
SO constants orders of magnitude higher than the ones 
obtained by first-principles calculations Q used here. 
In the experimental samples graphene was additionally 
coated by an AI2O3 layer to realize working tunnel bar- 
rier contacts. This likely brings metallic adatoms, which 
induce a stronger spin-orbit coupling strength, as has 
been reported for a full layer of Au atoms in contact 
with graphene, in which a several orders of magnitude 
larger BR-constant of about 13 meV was found [26[; sim- 
ilar large SO constants were predicted for impurities on 
graphene 27J. Suppose an adatom induces a local spin- 
orbit splitting of magnitude w 10 meV. The splitting 
spreads a distance s of perhaps a few bond lengths. Let 
the average distance between the randomly positioned 
adatoms be r. Then the DP spin relaxation rate is 
1/t s il 2 T(s/r) 2 . The rate is reduced from that for 
a homogeneous splitting by (s/r) 2 , which renormalizcs 
f2 2 due to the finite effective adatoms area. As a generic 
example we take s to be two bond lengths, s f=a 3 A, and 
a reasonable distance r « 10 nm, we get the spin relax- 
ation time t s ps 50 ps (using t ps 100 fs), being of the 
same order of magnitude as the measured value [2j. The 
adatom mechanism depends strongly on the adatom type 



and density, making it experimentally testable. 

In summary, we showed that charged impurities and 
polar optical surface phonons of the substrate gener- 
ate a random Bychkov-Rashba SO-field which leads to 
an almost temperature-independent spin relaxation in 
graphene. The calculated spin relaxation times give the 
upper bounds of what one can expect experimentally for 
a clean graphene on a substrate. The above calculation 
also shows that spin injection and spin transport should 
be severely limited if metallic electrodes are deposited 
directly on graphene. 

This work has been supported by the DFG SFB 689 
and SPP 1285. We thank E. Sherman for very valuable 
and inspiring discussions. 
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